arXiv: astro-ph/0112082v 1 4 Dec 2001 


Mon. Not. R. Astron. Soc. 000 , 000-000 (0000) Printed 1 February 2008 (MN style file vl.4) 


The dynamics of SO galaxies and their Tully-Fisher relation 


A. Mathieu , 1 M.R. Merrifield 1 and K. Kuijken 2 

1 School of Physics and Astronomy, University of Nottingham, Nottingham NG7 2RD 
2 Kapteyn Institute, Postbus 800, Groningen 9700 AV, the Netherlands 


Received date; accepted date 


ABSTRACT 

This paper investigates the detailed dynamical properties of a relatively homogeneous 
sample of disc-dominated SO galaxies, with a view to understanding their formation, 
evolution and structure. By using high signal-to-noise ratio long-slit spectra of edge-on 
systems, we have been able to reconstruct the complete line-of-sight velocity distribu¬ 
tions of stars along the galaxies’ major axes. From these data, we have derived both 
model distribution functions (the phase density of their stars) and the approximate 
form of their gravitational potentials. 

The derived distribution functions are all consistent with these galaxies being 
simple disc systems, with no evidence for a complex formation history. Essentially no 
correlation is found between the characteristic mass scale-lengths and the photometric 
scale-lengths in these galaxies, suggesting that they are dark-matter dominated even 
in their inner parts. Similarly, no correlation is found between the mass scale-lengths 
and asymptotic rotation speed, implying a wide range of dark matter halo properties. 

By comparing their asymptotic rotation speeds with their absolute magnitudes, we 
find that these SO galaxies are systematically offset from the Tully-Fisher relation for 
later-type galaxies. The offset in luminosity is what one would expect if star formation 
had been suddenly switched off a few Gyrs ago, consistent with a simple picture in 
which these SOs were created from ordinary later-type spirals which were stripped of 
their star-forming ISM when they encountered a dense cluster environment. 

Key words: Galaxies: lenticular - galaxies: structure - galaxies: kinematics and 
dynamics 


1 INTRODUCTION 

Since SO galaxies seem to have properties intermediate be¬ 
tween elliptical and spiral galaxies, Hubble (1936) placed 
these gas-poor disc systems between the spirals and ellipti¬ 
cals, at the bifurcation of his tuning fork galaxy classifica¬ 
tion scheme. However, three-quarters of a century later, it 
remains an open question as to whether this arrangement in 
any way reflects the physical origins of SO galaxies. 

One approach to addressing question of this kind has been to 
search directly for signs of evolution in galaxy morphologies 
by observing samples of these systems over a wide range of 
redshifts. In a classic study of this kind, Butcher & Oemler 
(1978) found that the fraction of blue galaxies in clusters 
was much higher in the past than it is now. They surmised 
that blue spirals were being converted to earlier-type sys¬ 
tems. This discovery fitted in extremely well with Gunn & 
Gott’s (1972) suggestion that SO galaxies could form from 
spiral galaxies in the dense environment of a cluster via tidal 
and ram pressure stripping. More recently, high-resolution 
observations (e.g. Couch et al. 1994; Lavery, Pierce & Mc¬ 
Clure 1992) have confirmed that the blue galaxies in high 


redshift clusters are, indeed, relatively normal spiral galax¬ 
ies, which could well be the progenitors of SO systems. 

The only problem with such analyses is that one is neces¬ 
sarily observing different galaxies in the nearby and distant 
samples, so any inference about the evolution from one type 
to another is of a circumstantial nature. A convincing case 
therefore requires that one look in some detail at the “fin¬ 
ished articles,” in order to find any archaeological evidence 
for the proposed evolution. Since galaxies are intrinsically 
dynamical entities, one should be able to extract impor¬ 
tant clues from their stellar kinematics, as derived from 
absorption-line spectra, as well as their photometric pro¬ 
perties. With high quality spectral data, one can measure 
the full line-of-sight velocity distribution (e.g. Koprolin & 
Zeilinger 2000, Fisher 1997), providing information not only 
on the motions of the stars, but also the gravitational po¬ 
tential responsible for the motions. 

For SO galaxies, one important piece of stellar-dynamical ev¬ 
idence would be provided by an analysis of the Tully-Fisher 
relation. In later-type disc galaxies, there is a strong correla¬ 
tion between circular rotation speed and optical luminosity 
(Tully & Fisher 1977), and this relationship becomes even 
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tighter when one looks in the near infrared (Pierce & Tully 
1992). If SO galaxies formed in a relatively benign way from 
spiral galaxies, one would expect that their optical luminosi¬ 
ties and circular rotation speeds would be little affected by 
the process, so a Tully-Fisher relation should still be appar¬ 
ent; the only significant difference would be that the strip¬ 
ping of their gas would switch off the star formation process, 
so the optical luminosities of the SOs should fade over time, 
shifting the zero-point of the relation. 

To-date, the most thorough search for a Tully-Fisher rela¬ 
tion in SO galaxies was that made by Neistein et al. (1999). 
This analysis is complicated by the fact that there is no sim¬ 
ple measure of the circular speed in an SO galaxy: there is 
no gas at large radii moving on circular orbits, and the stars 
follow significantly elliptical orbits, so their mean streaming 
velocity at any radius is lower than the local circular veloc¬ 
ity. Neistein et al. overcame the problem by appealing to 
the equations of galactic dynamics: by making a few sim¬ 
plifying assumptions, they were able to use the asymmetric 
drift equation (Binney & Tremaine 1987) to combine the 
mean streaming velocity and velocity dispersion of the stars 
in order to estimate the circular velocity at each radius. 
Although this analysis revealed some evidence for a trend 
between I-band luminosity and circular speed, Neisten et 
al. found a huge scatter in the relation, and that there was 
very little offset in the mean from the Tully-Fisher relation 
of later-type galaxies. They therefore concluded that these 
systems could not all have formed from the simple stripping 
of spiral galaxies. Instead, they suggested that the SO clas¬ 
sification actually represents a rather heterogeneous class 
of galaxies, which formed through a rather wide variety of 
processes. They further suggested that the absence of an off¬ 
set in the Tully-Fisher relation could be understood if the 
SO galaxies have more massive discs, so any fading in their 
luminosities is offset by the larger number of stars. 

In order to explore these conclusions a little further, this pa¬ 
per presents a detailed dynamical analysis of six SO galaxies. 
These galaxies have been selected to contain relatively small 
central bulges. If any SOs formed from simple stripping pro¬ 
cesses, one would expect their bulges to be little affected. 
Thus, by selecting SOs with the small bulges characteristic of 
later-type spirals, one might hope to pick out a relatively ho¬ 
mogeneous subsample of systems that formed via this route. 
Unfortunately, it is only when SO galaxies are very close to 
edge on that one can reliably determine that their bulges 
are small. As Neistein et al. (1999) demonstrated, the line- 
of-sight integration of starlight through such an edge-on disc 
means that there is quite a large correction to convert the ob¬ 
servable mean line-of-sight velocity into the circular stream¬ 
ing velocity of the stars. As a further complication, the ob¬ 
served velocity distribution for any line of sight through an 
edge-on disc will be highly non-Gaussian, due to the con¬ 
tribution from stars at large radii with small line-of-sight 
velocities. A dynamical analysis based on moments derived 
from a Gaussian fit to the line-of-sight velocity distribution 
is therefore prone to systematic error. To obviate these dif¬ 
ficulties, the current analysis uses data with a high enough 
signal-to-noise ratio for the complete line-of-sight velocity 
distribution to be derived, and these data are then fitted 
to a complete dynamical model in order to derive both the 
stellar distribution function and the gravitational potential 
needed for the Tully-Fisher relation. 


.♦ * 



Figure 1. Images from Digitized Sky Survey for the sample 
of SO galaxies. From Left to Right, Top to Bottom: NGC1184, 
NGC1611, NGC2612, NGC3896, NGC4179 and NGC5308. 

The remainder of the paper is laid out as follows. Section ^ 
describes the sample, the data analysis, the model fitting, 
and the resulting distribution functions and rotation curves. 
Section discusses the correlations in derived quantities, 
with particular reference to the Tully-Fisher relation. Sec¬ 
tion yl presents the conclusions that can be drawn from this 
analysis. 

2 DATA ANALYSIS 

2.1 The sample 

The galaxies chosen for this study were selected from the 
sample presented in Kuijken, Fisher & Merrifield (1996). 
This galaxy sample contains 28 early-type disc galaxies, and 
was originally observed to search for counter-rotating popu¬ 
lations, so high signal-to-noise ratio spectra were obtained. 
Details of the spectroscopic observations and standard data 
reduction can be found in Kuijken, Fisher & Merrifield 
(1996). From this sample we selected the six SO galaxies 
of with similar mean rotation speeds, so that there is some 
chance that the systems are comparable in their dynamics. 
We also selected galaxies that are close enough to edge-on 
for it to be apparent that the systems are disc-dominated, 
with relatively small bulges (and thus that a simple evo¬ 
lutionary path from later-type spiral galaxies is plausible). 
Images of the sample galaxies are shown in Fig. 111. 

2.2 Kinematic analysis 

LOSVDs were extracted from the reduced spectra using the 
unresolved gaussian decomposition (UGD) algorithm (Kui¬ 
jken & Merrifield 1993). This method models each LOSVD 
as the sum of a set of unresolved gaussians with fixed means 
and dispersions. The amplitudes of the gaussians are varied 
so as to minimize the difference between the galaxy spec¬ 
trum and the model derived by convolving the LOSVD with 
a suitable template star spectrum. This approach has the 
advantage that it does not force any particular functional 
form on the LOSVD, beyond the smoothness imposed by 
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Figure 2. Top left: The line-of-sight velocity distribution as a 
function of projected radius along the major axis of NGC 1184. 
Top right: the corresponding quantities for the best-fit dynamical 
model. Bottom left: plot of [(data-model)/<r] 2 , showing where the 
major y 2 residuals arise. Bottom right: ln(y 2 ) as a function of 
the adopted parameters for the gravitational potential, (ro,ro); 
the solid contour shows where ln(y 2 ) = ln(l.l X Xmin )■ 


the widths of the component gaussians. It is therefore well- 
suited to modeling complex systems like edge-on disc galax¬ 
ies, where integration along the line of sight, potentially 
through multiple components, can make the shape of the 
LOSVD rather complex. 

In order to increase the signal-to-noise ratio of the resulting 
kinematic data, the two-dimensional LOSVD as a function 
of position along the major axis, F(vi as , R), was assumed to 
have the symmetry of an edge-on axisymmetric disc system, 
so that F(v ios, R) = F(—vi os , —R). We therefore determined 
the kinematic coordinates of the centre of each galaxy by 
finding the point in F(vi oa ,R) at which this symmetry is 
most closely obeyed. We then averaged the data from the 
two sides of each galaxy. The resulting mean estimates for 
F(v\ oa ,R) are shown in the top left panels of Figs. 110 
0 | and 0 


2.3 Dynamical analysis 

Having derived these kinematic data, we must now attempt 
to describe them using a dynamical model. The simplest 
credible model for the major axis kinematics of a galaxy in¬ 
volved fitting them to a two-integral distribution function, 
f(E,L), which describes the phase-space density of stars, 
where E is the energy and L is the angular momentum of the 
stars about the symmetry axis, two of the integrals of motion 
(Binney & Tremaine 1987). Such a model would be exact for 
an infinitely-thin axisymmetric disc, but it is also a credi¬ 
ble approximation for a system of finite thickness. Defining 
the usual polar coordinates, the energy in the z direction, 
E z = &(R, z ) — <&(72, 0) + |u 2 , is approximately an integral 


Figure 3. Same as in Fig. 


for NGC 1611 



Figure 4. Same as in Fig. 


for NGC 2612 


of motion for a thin disc. A system whose distribution func¬ 
tion takes the form f(E, L)g(E z ) will have dynamics in the 
z — 0 plane identical to those of an infinitely-thin disc with 
a distribution function f(E,L). Thus, the dynamics in the 
plane of the disc can credibly be modeled by treating the 
major-axis kinematics as if they were those of an infinitely- 
thin disc. 

Mathieu & Merrifield (2000) presented an algorithm by 
which the observed density of stars as a function of pro¬ 
jected radius and line-of-sight velocity, F(R,v i os ) could 
be iteratively inverted to find the distribution function, 
f(E,L), that would produce such observable kinematics. 
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Figure 5. Same as in Fig. 


for NGC 3986 
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Figure 7. Same as in Fig. 


for NGC 5308 



Figure 6. Same as in Fig. 


for NGC 4179 


This method builds an estimate of the distribution func¬ 
tion using just the data from the upper part of the line-of- 
sight velocity distribution where the velocities exceed the 
circular velocity. If the correct gravitational potential is 
adopted, then the rest of the velocity distribution automat¬ 
ically matches the data; however, a mismatch will occur if 
the wrong potential is assumed. Thus, this algorithm algo¬ 
rithm returns not only an estimate for the disc distribution 
function, but also constrains the form of the gravitational 
potential. 

In principle, the gravitational potential in the disc plane 
could take any form, but the data are not of high enough 


Name 

vo [km/s] 

ro [arcsec] 

rd [arcsec ] 

Mj 

M h 

NGC 1184 

271 ± 17 

1.4 ±0.4 

21.3 ± 1.5 

-21.1 

-23.0 

NGC 1611 

271 ± 15 

1.8 ±0.3 

12.1 ±0.3 

-21.5 

-23.4 

NGC 2612 

252 ± 11 

2.2 ±0.3 

12.8 ±0.2 

-19.9 

-21.8 

NGC 3986 

289 ± 06 

2.4 ±0.6 

07.6 ± 1.1 

-21.6 

-23.5 

NGC 4179 

275 ± 07 

1.7 ±0.4 

12.8 ±0.5 

-20.9 

-22.8 

NGC 5308 

321 ± 17 

0.5 ±0.1 

16.5 ±0.4 

-21.6 

-23.5 


Table 1. Estimates of the parameters (ro,no) for the gravita¬ 
tional potential, estimate of the photometric disc scale length r,/ 
and absolute I and H magnitudes for each galaxy. 


quality to allow a completely non-parametric derivation. We 
therefore adopt the simple but sufficiently general form of a 
softened isothermal sphere potential of the form 




(1) 


For this potential, the circular velocity can be written 


v c {r) 


vp r 

\A 2 + r o ’ 


( 2 ) 


with the parameters ro and vo specifying how quickly the 
rotation curve rises and its asymptotic value. 

Using the iterative algorithm, distribution functions have 
been derived for all the sample galaxies. Figure ^ shows 
the best-fit distribution functions obtained in this way. All 
the systems show the single concentration of low angular- 
momentum material, characteristic of a disc. Thus, there is 
no evidence even in these detailed dynamical models of any 
peculiarities in the stellar dynamics that one might expect 
if these systems had formed in a spectacular manner such 
as in a merger (Bekki 1998). 

To assess the reliability of the fitting process, Figs. [| [I], 
|E| |g| and |ij show the “observable” kinematics of these dy¬ 
namical model, F(R,v i os ), and the difference between these 
predictions and the observed kinematics (normalized by the 
uncertainty in the observations). Finally, these figures also 
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Figure 8. The disc distribution functions for the galaxies as derived iteratively from their line-of-sight kinematics using the best-fit 
gravitational potentials. 



Figure 9. Same as in Fig. for NGC 3986, using data on the Figure 10. Same as in Fig. for NGC 3986, using data on the 

’’left” side of the galaxy. ’’right” side of the galaxy. 


show how rapidly the goodness of fit (in a x 2 sense) degrades 
as one goes away from the optimal values for ro and Vo in 
the gravitational potential. The best-fit values for ro and Vo 
are listed in Table [I| 

So far in our analysis we assumed that the galaxies are per¬ 
fectly axisymmetric and we used average data from both 
sides of the galaxy. However the kinematics of these galaxies 
show small asymmetries and we also performed the fitting 
process on each side of the galaxies separately. We use the 
best-fit values of r o and Vo on each side of the galaxy to 


estimate the dispersion of r o and Vo around the mean for 
each galaxy. An example of fits on both sides of one galaxy 
(NGC3986) is shown in Figs. and [u| 

Clearly, in these galaxies there are some systematic residuals 
in the difference between model and data. Indeed, the mini¬ 
mum values of x 2 are well above what one might expect for a 
formally good fit. However, given the assumptions involved 
as to the separable form of the distribution function and the 
adopted parametric function for the gravitational potential, 
these small residuals are not particularly surprising. In or- 
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Figure 11. Same as in Fig. y for a bulge-disc model of a 
NGC 1184 


der to assess their impact on our subsequent analysis, we 
have explored more sophisticated models for the distribution 
function. In particular, we have investigated what happens 
if our assumption of a single-component disc-like distribu¬ 
tion function is relaxed, by explicitly adding the contribution 
from a central bulge component (rather than treating it as 
part of the disc). To this end, we adopt a simple bulge model 
with a distribution function of the form 

f(E, L ) = a[exp(-(3E) - exp(-/J£o)] exp( 7 L), (3) 

where (E < Eg) is the energy of a star and L its total 
angular momentum, and where a, /3 and 7 are constants 
(see Jarvis & Freeman, 1985). We then computed the line- 
of-sight velocity distribution corresponding to this distribu¬ 
tion function and subtract it from the observed line-of-sight 
velocity distribution. The choice of the constants remains 
somewhat arbitrary as a full bulge-disc decomposition is 
underconstrained, but we stick to values of the parameters 
such that the estimate of line-of-sight velocity distribution 
of the disc alone (after subtracting the line-of-sight veloc¬ 
ity distribution produced by the bulge from the observed 
line-of-sight velocity distribution) is positive almost every¬ 
where. Rather surprisingly, this extra spherical component 
makes very little difference to the best-fit gravitational po¬ 
tential: Fig. [ill illustrates the best-fit model for NGC 1184, 
in which the constraints on ro and vo are close to those in 
the pure disc model. It would thus appear that the simple 
model adopted above does a robust job of estimating the 
gravitational potential even where the assumption of a thin 
disc is not strictly valid. 


2.4 Photometric analysis 

B-band integrated magnitudes for the sample galaxies have 
been taken from the RC3 catalogue, where available, or from 
the Lyon-Meudon extragalactic database (LEDA) otherwise. 


These have been converted to absolute magnitudes using dis¬ 
tance estimates from the LEDA database, where available. 
For NGC1611 and NGC2612, as no distance is available, 
we estimate the distance using the systemic velocity derived 
from our spectra assuming Ho = 80 kms~ x Mpc -1 . It would 
be useful to study the integrated photometric properties of 
these systems in a number of bands. Unfortunately, such 
data are unavailable for most systems. However, the spread 
in the integrated colors of SO galaxies is very small: for the 
SO galaxies in the sample of Persson, Frogel & Aaronson 
(1979), B — H has a mean value of 3.9 and a scatter of only 
~ 0.3, while B — I has a mean value of 2 and a scatter of 
~ 0.3 for the SO galaxies in the sample of Neistein et al. 
(1999). As we shall see below, these scatters are comparable 
to other sources of uncertainty, so for this analysis we can 
simply apply the mean color corrections in order to estimate 
absolute magnitudes in the I and H bands. The resulting 
magnitudes are listed in Table l]. 

The principal spatially-resolved photometric quantity of in¬ 
terest in these disc-dominated systems is the disc scale- 
length. To determine this quantity, we first estimate the in¬ 
tensity profiles, 7(r), as a function of radius by integrating 
the line-of-sight velocity distributions derived above with re¬ 
spect to velocity. We adopt this approach in preference to 
using imaging data, as it guarantees that the photometric 
and kinematic measurements are probing the same stellar 
population. We then approximate the starlight of the disc 
at large radius by fitting a function of the form 

Jmod(r) = 2Ad\r\K\(\r\/rd), (4) 

where Ad and rd are free parameters, and A'i is a modified 
Bessel function, which represents the projected density pro¬ 
file of an edge-on exponential disc. With only rather noisy 
photometric profiles and no detailed surface photometry for 
our galaxies, we only fit a disc model at large radius as fit¬ 
ting a bulge-disc model proves very much underconstrained 
and does not give reliable values for the parameters of the 
disc and bulge. The best fit photometric models are shown 
in Fig. [Ijj where the models are fitted independently on 
each side of the galaxy, and the derived photometric scale- 
lengths rd are given in Table |I|. This quantity differs from 
the photometric scale length usually estimated by fitting 
two-component bulge-disc models . 


3 DISCUSSION 

Having derived the dynamical and photometric properties of 
these galaxies, we are now in a position to look for clues to 
their origins. Figure [jj| shows the physical parameters of the 
rotation curves in these galaxies. There are no signs of any 
correlation between ro and Vo, suggesting that these galaxies 
occupy dark halos with the usual wide range of characteris¬ 
tics, although the sample is really too small to say anything 
definitive. Similarly, there is no discernable correlation be¬ 
tween the mass scale-length, r 0 , and the photometric scale- 
length, rd, (see Fig. [ll]) arguing against Neistein et al .’s 
(1999) suggestion that these galaxies should be more dom¬ 
inated by the mass in their stellar discs than normal spiral 
galaxies. 

A more interesting result comes when we look at the Tully- 
Fisher relation for SO galaxies. Figure [1^ shows the I-band 
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Figure 12. Plot of the logarithm of the projected luminosity density as a function of projected radius along the major axes of the six 
sample galaxies, as derived from their line-of-sight velocity distributions. The overplotted smooth lines show the best-fit disc model. 
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Figure 13. Plot of the estimates of the maximum circular veloc¬ 
ity vo vs. the estimates of the parameter ro- 


relation derived by Neistein et al. (1999), which lie quite 
close to the Tully-Fisher relation for normal spiral galaxies 
(Pierce & Tully 1992), but with a large RMS scatter of ~ 0.6 
magnitudes about the best-fit line parallel to the standard 
Tully-Fisher relation. The galaxies in the current sample lie 
systematically below the Neistein et al. galaxies. With the 
galaxies all selected to have comparable rotation speeds, we 
clearly cannot derive a slope for the Tully-Fisher relation 
from these data, but fixing the slope to that of the relation 
for normal spiral galaxies, we find that the best-fit line is off¬ 
set from the spiral galaxy relation by ~ 1.8 magnitude, with 
a scatter about the best-fit line of only ~ 0.3 magnitudes. 
For this fit we exclude NGC2612 as its distance estimate is 
rather uncertain due to a small systemic velocity. Even al¬ 
lowing for the small sample size, this scatter is smaller than 



0 5 10 15 20 25 

r d [arcsec| 

Figure 14. Plot of the estimates of photometric disc scale length 
r,i vs. the estimates of the parameter ro for the gravitational 
potential. 


that in the Neistein et al. sample at a statistically significant 
level. 

This small scatter and offset in the Tully-Fisher relation by 
~ 1.8 magnitude in the /-band is what one would expect if 
star formation had been suddenly switched off a few Gyrs 
ago so that these SO galaxies contained just the old stellar 
population of a normal spiral galaxy. Using stellar popula¬ 
tion models (Chariot & Bruzual 1991), assuming SOs galax¬ 
ies had had similar formation histories as late-type spirals 
until a few Gyrs ago, we would expect the stellar population 
of SOs to have faded significantly, resulting in a noticeable 
offset in I magnitude, such as the one observed in our sam¬ 
ple. 

This point is made even more dramatically if we convert 
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Figure 15. Absolute /-band magnitude versus log^o/km/s) for 
the sample of SO galaxies in Neistein et al. (1999), and those in 
this analysis (filled triangles). The line shows the /-band Tully- 
Fisher relation for later-type spirals from Pierce &; Tully (1992). 
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Figure 16. Estimated absolute l /-band magnitude versus 
log (no/km / s) for the sample of SO galaxies in Neistein et al. 
(1999), and those in this analysis (filled triangles). The line shows 
the //-band Tully-Fisher relation for later-type spirals from Pierce 
& Tully (1992). 

to estimated IT-band magnitudes using the prescription in 
Section [^. In this band, the luminosity is dominated by the 
old stellar populations; as Fig. [l6] shows, the SOs in the cur¬ 
rent sample lie quite close to the Tully-Fisher relation for 
later-type spiral galaxies, suggesting that their old stellar 
populations are rather similar. 


4 CONCLUSIONS 

In this paper, we have calculated the first detailed dynam¬ 
ical models for a small sample of edge-on SO galaxies with 
small bulges. In addition to producing distribution functions 
for these disc-dominated systems, which look very much as 
one would expect for normal disc systems, the analysis also 
returned estimates for the parameters of their gravitational 
potentials. 

The interpretation of these data all points to a simple pic¬ 
ture in which these systems were formed by the stripping of 


gas from normal spiral galaxies. The distribution functions 
are all well modeled by unexceptional stellar discs, similar 
to those expected in the old stellar populations of spiral 
galaxies. In addition, the galaxies obey a reasonably tight 
Tully-Fisher relation, which is offset from the relation for 
normal spiral galaxies by the amount that one would expect 
if star formation had been shut off a few Gyrs ago, so that 
all that remains in these systems are the rather fainter old 
stellar populations. 

This result appears to conflict with Neistein et al.’s (1999) 
analysis, which showed a much greater scatter in the Tully- 
Fisher relation with less systematic offset. Part of the differ¬ 
ence may be due to the lower signal-to-noise ratio of the data 
in their larger sample, which limited their ability to carry out 
detailed dynamical modeling, particularly for galaxies that 
lie very close to edge-on. However, there is also a systematic 
difference in the way that the samples were selected: the 
edge-on galaxies in the analysis of this paper were specif¬ 
ically chosen to contain small bulges. This selection crite¬ 
rion means that these galaxies are prime candidates to have 
formed from gas-stripped spiral galaxies. If, as Neistein et 
al. suggest, SOs are a “mixed bag” that formed in a variety 
of ways, it should come as no surprise that this particular 
subsample obey a tight Tully-Fisher relation that is not seen 
in the general population of SOs. To test such hypotheses, 
we ultimately need data of the quality presented in this pa¬ 
per for a much larger sample of galaxies, extending both the 
range of absolute magnitudes and of other parameters such 
as the bulge size. 
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